Disorder-induced pseudodiffusive transport in graphene nanoribbons 
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We study the transition from ballistic to diffusive and localized transport in graphene nanoribbons 
in the presence of binary disorder, which can be generated by chemical adsorbates or substitutional 
doping. We show that the interplay between the induced average doping (arising from the non-zero 
average of the disorder) and impurity scattering modifies the traditional picture of phase-coherent 
transport. Close to the Dirac point, intrinsic evanescent modes produced by the impurities dominate 
transport at short lengths, giving rise to a new regime analogous to pseudodiffusive transport in clean 
graphene, but without the requirement of heavily doped contacts. This intrinsic pseudodiffusive 
regime precedes the traditional ballistic, diffusive and localized regimes. The last two regimes 
exhibit a strongly modified effective number of propagating modes, and a mean free path which 
becomes anomalously large close to the Dirac point. 



I. INTRODUCTION 

Graphene continues to fascinate, due in particular to 
its transport properties close to the Dirac point, where 
a finite conductivity is observed experimentallyi" in spite 
of the vanishing density of states. For clean graphene, 
simple Dirac fermion models^ prove to be sufficient 
to describe this effect based on a gapless spectrum of 
evanescent modes, pinned at the reservoir contacts, that 
allow for quantum tunneling through macroscopically 
large graphene samples! 3 ^ Surprisingly, the full transport 
statistics becomes indistinguishable from that of diffu- 
sive metalsj^for which reason the name 'pseudodiffusive' 
was coined to describe this transport regime. Further- 
more, these evanescent modes can give rise to Hanbury 
Brown-Twiss type cross-correlations in a ballistic multi- 
probe graphene structure,^ similar to those observed in 
conventional diffusive conductors. 

Just as ballistic transport, the effects of disorder in 
graphene display a surprising richness. E.g., armchair 
ribbons are found to be more sensitive to bulk disorder as 
their zigzag counterparts, while edge disorder can open 
up transport gaps in both of themP In general, differ- 
ent types of disorder can be classified in terms of pre- 
served symmetries, including chirality which is preserved 
for disor der tha t is smooth on the scale of interatomic 
distance d 7 * 8 * 9 * 10 ! and significantly modifies the standard 
results for transport in disordered metals. For example, 
carriers in graphene ribbons with smooth disorder cannot 
be localized,^ and their conductivity is seen to increase 
with system size L as a ~ y/L/ l m in ID j^j an d ~ hxL/l m 
in 2D, where l m is the mean free pathP^H] Short-range 
disorder, on the other hand, does not preserve chiral- 
ity, and is therefore generally understood to induce a 
sequence of ballistic transport, diffusion, and An derson 
localization in quasi-lD wires of increasing lengthP 15 * 16 * 

Until now, the influence of evanescent modes close to 
the Dirac point has been studied in ballistic systems only, 
and it is therefore unclear whether their peculiar spec- 



W 




L 



FIG. 1: (Color online) Schematic view of an armchair ribbon 
with width W . Disorder is introduced in a region of length 
L by adding a potential Vb to a percentage p of lattice sites. 
Disorder sites are depicted by red dots. The induced local 
doping, due to the non-zero average of the disorder, is effec- 
tively described by a barrier of height V = p Vb . 



trum could manifest itself in some significant way in the 
conductivity of disordered graphene nanoribbons. In this 
work we will show that evanescent modes play a cru- 
cial role for the transport in graphene nanoribbons with 
short-range disorder. In particular, intrinsic evanescent 
modes arising from the disorder itself (i.e., independently 
of the contacts) can reinstate an 'intrinsic pseudodiffu- 
sive' regime which precedes the ballistic, diffusive and 
localized regimes. This new regime arises if the short- 
range disorder does not vanish on average and therefore 
induces systematic doping, as is the case, e^g., for chemi- 
cal adsorbates^ or substitutional doping We describe 
such systems using a binary disorder model (the random 
binary allo y 19 ! 20 ') and obtain our results by combining 
numerical computations with analytical considerations. 
For longer ribbons, the intrinsic pseudodiffusive regime 
crosses over into the traditional diffusive and localized 
regimes, but we find that the mean free path becomes 
anomalously large at an effective Dirac energy which 
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takes the induced doping into account. 



II. BINARY DISORDER MODEL 

We base our investigation of binary disorder in a 
graphene nanoribbon on the following tight-binding 
Hamiltonian, 
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The first term describes clean graphene with hopping 
amplitude t rs 2.9eV between nearest neighbors on the 
hexagonal lattice of carbon atoms (with interatomic dis- 
tance a = 1.49 A). The second term describes binary 
disorder, generated by onsite energies Vi which are un- 
corrected random variables that take values i>hi g h = Vq 
with probability p, and i>i ow = with probability 1 — p. 
We have dropped spin indices in the model since spin is 
a passive degree of freedom in all our discussion. 

Physically, this model describes, for example, a 
graphene nanoribbon with a fraction p of its sites cou- 
pled to identical chemical species adsorbed on its sur- 
face. If these adsorbates are screened their presence is 
equivalent to a (real and energy dependent) self-energy 
term acting on each of the carbon sites coupled to an 
adsorbateP^ Typical values of Vq can be estimated from 
DFT calculations,^ and for species such as Ag, Cu or 
Au, one finds values Vo » 0.5eV « 0.2t which are small 
as compared to t. 



III. RESULTS AND DISCUSSION 

To probe the intrinsic transport properties in the pres- 
ence of binary disorder we present results of numerical 
computations on nanoribbons of width W and a dis- 
ordered region of length L, with the reservoir contacts 
placed at x = ±oo (see Fig. [I]). These computati ons are 
based on the recursive Green's function technique ! 22 * 23 * 24 ^ 
We obtain qualitatively similar results for metallic and 
semiconducting armchair ribbons as well as for zigzag 
ribbons; we therefore only show results for metallic arm- 
chair ribbons. 

Figure [2] shows the averaged conductance of a ribbon 
of width W = 103^30/2 as a function of the length of the 
disordered region; different curves correspond to different 
values of the chemical potential. The disorder parameters 
used, p = 0.4 and Vq — 0.25i, correspond to a substantial 
concentration of adsorbates that are weakly hybridized 
with their carbon hosts; a sample size of 50 impurity 
configurations is used in the average. 

Two distinct transport regimes immediately catch the 
eye. The conductance first decreases sharply over a 
length scale comparable to the ribbon width W, and then 
abruptly crosses over to a much slower decay. The con- 
ductivity a = (G)L/W in the pseudodiffusive regime, 
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FIG. 2: (Color online) Averaged conductance (G) (in units 

Go = ^f-) of an armchair ribbon of width W = 103\/3ffl/2 as 
a function of the length of the disordered region for different 
chemical potentials. The parameters of the binary disorder 
are p = 0.4, Vo = 0.25t, so that the effective Dirac point in 
the disordered region is shifted from zero to V = pVo = O.lt. 
Dashed curves are fits to Eq. j3j, with l m as a fit parame- 
ter. Inset: Pseudodiffusive regime zoom-in. Thin grey curves 
correspond to different disorder realizations. Good agreement 
is obtained with the ballistic conductance through a square 
potential barrier of height V, shown in dashed black, see Eq. 



shown in Fig. [3] develops a clear plateau at the value a = 
2Gq/tt (the minimal conductivity in clean graphene^') 
when the chemical potential fi coincides with the average 
value V = pVo of the disorder. We now argue that this 
transport regime indeed constitutes an intrinsic pseudod- 
iffusive regime arising from disorder-induced evanescent 
modes, while the crossover to the truly diffusive regime is 
characterized by a mean free path which becomes anoma- 
lously large at = V. 



A. Intrinsic pseudodiffusive regime 

To understand the numerical results for short lengths 
it is useful to compare the average conductance of Fig. 
|2] to the conductance of the same nanoribbon in which 
the disordered region is replaced by a potential barrier 
of height equal to the average disorder doping V — pVq. 
We compute the conductance through the barrier, Gb, by 
wave-matching within the Dirac equation description of 
low-energy transport, which yields 



Gb — Go 



E 



(1 - z 1lk f (1 - z nk >) 2 



(znk 



\ 2 i —ik'L ( 1 \ 



(2) 

Here the conductance quantum is defined by Go = 2e 2 //i, 
the sum runs from n = -(N - l)/2 to (N - l)/2, 
the number of incoming propagating modes is N — 
1 + 2Int[fj,W/(Trhv)], fi is the reservoir chemical po- 
tential, the Fermi velocity is given by v — ^ta/h, 
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(k + iq n )/y/k 2 + ql, k = yj {fi/hv) - q%, k' 

\ [(jU - V)/hv] 2 - ql, and q n = mr/W. 

The comparison between G and G& over the range of 
the abrupt decay is shown in the inset of Fig. [2] where 
each grey curve corresponds to an individual disorder 
configuration, and the dashed black curve corresponds to 
Eq. ([2]) . We obtain a remarkably good agreement of both 
curves, even down to hallmark features such as Fabry- 
Perot oscillations. This good agreement only begins to 
break down at lengths for which the disorder starts to 
mix modes, i.e., the mean free path l m , which will be 
discussed below. 

The implications are clear: The sharp decline of con- 
ductance observed at lengths L < W is due to the decay 
of evanescent modes induced by the sudden change in 
local average potential V. This average potential shift 
is well defined because the effect of the disorder is self- 
averaging at the considered impurity concentrations.. Al- 
though these evanescent waves are produced by the im- 
purities, their combined effect is equivalent to evanescent 
modes pinned at the boundary of the disorder region, 
x = and x = L, which allow tunneling through an 
effective potential barrier V = pVq. Furthermore, their 
contribution to transport at energies fi ~ V and very 
wide ribbons produces the universal statistics of diffu- 
sive metals. Therefore, the observed fast decay in aver- 
age conductance is the signature of a disorder-induced 
pseudodiffusive regime. As we will show in the following 
section, the self averaging of the disorder arises since the 
mean free path becomes anomalously large around the 
effective Dirac point /u = V. Contrary to the ballistic 
regime, the intrinsic pseudodiffusive regime is therefore 
universal in the limit of wide and short ribbons, which is 
verified by the negligible spread of the curves in the inset 
of Fig. [2] over the range of the abrupt decay. 



B. Ballistic to diffusive crossover 

At lengths L ~ W evanescent modes through the dis- 
ordered region have decayed, and give way to ballistic 
propagation of the remaining N p (fi) = 1 + 2Int[(/i — 
V)W/(nhv)] modes, up to lengths L comparable to the 
mean free path l m , where the traditional diffusive regime 
sets in. In the range W < L < N p l m , therefore, the 
large- N p form of the average conductance in the ballistic 
to diffusive crossover becomes approximately validj^ 
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This is followed by a transition into a localized regime for 
nanoribbon lengths greater that the localization length 

decays exponentially as exp(— 2L/£) 

The mean free path in the bulk graphene sample is 
known to be inversely proportional to the electron energy 



N p l m , for which (G) < Go- From that point on, (G) 
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FIG. 3: (Color online) Conductivity a = {G)L/W versus 
length for an armchair ribbon of width W = 103\/Sa/2 at 
the effective Dirac point fi = V and fixed impurity fraction 
p = 0.4. Different curves correspond to different strength 
Vo — V/p of the disorder. As the chemical potential of the 
effective Dirac point increases, the number N of incoming 
modes grows, and a plateau develops at the universal conduc- 
tivity minimum 2Go/^ (dashed black line). 
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FIG. 4: (Color online) Mean free path as a function of chem- 
ical potential for a metallic armchair nanoribbon of width 
W — 103\/3a./2, impurity strength Vo = 0.25£, and impurity 
fraction p — 0.4. Black curve: bulk value, Eq. Q. Red curve: 
finite width result, Eq. J5I. Dots: values obtained by fitting 
numerical results to Eq. (pi. 



and the disorder strength, see e.g. Refs. I26I27I Adapt- 
ing the results of these references to our binary disorder 
model, we arrive at the following expression for the mean 
free path 



3a/3 at 3 
4 p(l-p)V*\n- 
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The — V\ divergence can be traced back to the van- 
ishing bulk density of states of graphene at /i ~ V in 
the presence of average doping V = pVq. This depen- 
dence is consistent with the behavior observed in Fig. [2] 
for lengths greater than W. However, a detailed analysis 
shows that Eq. Q provides only a qualitative descrip- 
tion of our data. In order to improve the agreement one 
has to account for the finite size effects in the density of 
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FIG. 5: (Color online) Conductance of a graphene armchair 
ribbon of length L = 3000a and width W = 103^^/2 for 
different values of the impurity fraction p and scattering po- 
tential Vq. The chemical potential is fixed to fi = O.li. The 
solid line corresponds to the condition p, = V = pVo- The 
dotted lines are guides for the eye, bounding the region of 
increased conductance around the effective Dirac energy. 



calculated the conductance of a long metallic graphene 
ribbon (i = 3000a) for different parameters p and Vq, 
at a fixed chemical potential p = O.li. The results are 
summarized in Fig. [5] Parameters p and Vq that lead to a 
large mean free path show up as regions with a relatively 
high conductance. E.g., in the region around pVo = p 
depicted with a solid black curve in Fig. [5] the conduc- 
tance is close to G = 2e 2 /h, even for this rather long 
ribbon. The conductance is also large around p = 
and p = 1, which correspond to the regimes where the 
effective disorder strength p(l — p)V Q 2 is small. This fac- 
tor thus accounts for the semilunar shape of the phase 
diagram. Nevertheless, the physical picture leading to 
a high conductance is somewhat different for small and 
high disorder percentages. For small p (and Vo), we are 
in the strict ballistic regime, where the electrons do not 
feel the impurities at all, not even their average doping 
effect. On the other hand, for high p, we have effectively 
a disorder-induced potential barrier in the system. 



states. In the simplest weak disorder approximation the 
mean free path acquires the forrrP^ 
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where the effective density of states p(fi) is defined as 
follows 



pQi) = 1 + 2H£ 



(6) 



Here E n = hvnn/W are threshold energies at which new 
channels open up. According to Eq. ^ the mean free 
path is no longer divergent at fi = V, rather it saturates 
at a finite value Z max oc W. In Fig. 3] we show the com- 
parison of the mean free path Eq. (15 1 to the one resulting 
from fitting our data in Fig. [2] with Eq. we find that 
Eq. (|5| reproduces our numerical data very accurately. 



C. Dependence of conductance on disorder 
parameters 

In order to give a comprehensive overview of the inter- 
play between the various disorder parameters, we have 



IV. CONCLUSION 

In conclusion, we considered the conductance of 
graphene nanoribbons containing binary disorder (typ- 
ical for chemical doping) which induces a shift V of the 
effective Dirac point, along with scattering on single dis- 
crete impurities. For short ribbons we identified a trans- 
port regime dominated by evanescent modes which give 
rise to an intrinsic form of pseudodiffusive transport. At 
a longer length scale, a ballistic, and then a diffusive 
regime is entered where the mean free path is determined 
both by scattering on the individual impurities, as well 
as by the shifted potential. Close to the effective Dirac 
point, the mean free path becomes anomalously large, 
which ensures the universality of the intrinsic pseudodif- 
fusive regime. Compared to ordinary disordered conduc- 
tors, the presence of non-zero-average short-range disor- 
der typical for chemically doped graphene therefore adds 
a new stage to the conventional path towards localization. 
The predicted intrinsic pseudodiffusion should be di- 
rectly measurable in chemically functionalizcd graphene 
nanoribbons. 
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